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Abstract 

The compressible Navier-Stokes-Poisson system is concerned in the present paper, 
and the global existence and uniqueness of the strong solution is shown in the 
framework of hybrid Besov spaces in three and higher dimensions. 

Key words: compressible Navier-Stokes-Poisson equations, global existence and 
uniqueness, hybrid Besov spaces 



1 Introduction 



In the present paper, we consider the Cauchy problem of the following 
compressible Navier-Stokes-Poisson equations 

Pt + div(pu) = 0, 

(pu)j + div(pu (g) u) + VP(p) = pV0 + pAu + (p + A) V divu, 
A0 = p-p, 
. (p,u)(0) = (po,uo), 

for (t, x) G [0, +oo) X M^, ^ 3. p, u and (p denote the electron density, 
electron velocity and the electrostatic potential, respectively. -P(p) = is 
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the pressure with 7 = 2. A are the constant viscosity coefficients satisfying 
/i > and 2fi + A^A ^ 0. The constant p stands for the density of positively 
charged background ions. The Navier-Stokes-Poisson system is a simphfied 
model (for instance, the energy equation is not taken into granted) to de- 
scribe the dynamics of a charge transport where the compressible charged 
fluid interacts with its own electric field against a charged ion background [6]. 

Recently, many interesting researches have been devoted to many topics 
of the compressible Navier-Stokes-Poisson (NSP) system. The global existence 
of weak solutions of the compressible NSP system subject to large initial data 
is shown [8,19]. The quasi-neutral limits and related combining asymptotical 
limits are proven [7,9,13,18]. In the case that the potential force representing 
the self-gravity in stellar gases, the global existence of weak solutions and 
asymptotical behaviors are also investigated recently, and the stability analysis 
for compressible Navier-Stokes-Poisson and related systems is also carried out, 
refer for instance to [10-12,15] and references therein. 

The global existence of the classical solution is shown recently [14] in 
terms of the framework by Matsumura-Nishida. In addition, the influence of 
the electric field is justified, which affects the dissipation of the viscosity and 
the time-decay rate of global solutions of IVP (1.1) to the equilibrium state 
(p,0), namely, 

ci(l + t)-t < ||(p-p)(t)|U2(K3) <C(l+t)-t, (1.2) 
ci{l + ty-* < ||m(t)||i2(K3) < C(l + t)-3, (1.3) 

where the decay rate of the momentum or the velocity is slower than the rate 

3 

(l-|-t)~4 for compressible Navier-Stokes equations. A natural question follows 
then, that is, whether the similar phenomena can be shown for global weak 
solutions or strong solutions with lower regularity. 

To this end, the first step is to show the global existence of strong solutions 
in some Besov space with lower regularity. In this paper, with the help of the 
classical Friedrichs' regularization method, Littlewood-Paley analysis and hy- 
brid Besov spaces, we are able to construct the approximate solutions, obtain 
the a-priori estimates in hybrid Besov spaces, and prove the global existence 
of the unique strong solution by the compactness arguments as in [2,4,16]. 
Indeed, in terms of the div-curl decomposition we can decompose the velocity 
vector field into a vector field of the compressible part and a incompressible 
part. Then, the original compressible system for the density and the velocity 
can be decoupled into a system involving only the compressible system for the 
ir-rotational (compressible) part of velocity vector field and the electron den- 
sity and the diffusion equation for the divergence free (incompressible) part of 
velocity vector field as used in [4]. Thus, we can investigate the compressible 
velocity field part and the incompressible velocity field part separately to get 
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the expected estimates in some hybrid Besov spaces. As one can see later, 
however, the appearance of the electric field leads to the rotational coupling 
effect and the loss of regularity of density and the velocity vector field. 

For simplicity, we only deal with the case 7 = 2, the arguments used here 
can be applied to show the global existence for general 7 > 1. We have the 
main theorem as follows. 

Theorem 1.1 Let A^^3, 7 = 2, yu>0 and 2/i + A^A ^ 0. Assume po — p E 

~ K — l E. ~ £L — i £L — i 

^2,1 ^ ' ^ 0''nd uo £ -82^1 ^ ' ^ • Then, there exist two positive constants a small 
enough and M such that if 

\\Pn — P\\ N 5 JV + Un JV 3 JV , ^ Q!, 
^2,1 ^2,1 

then (1.1) yields a unique global solution (p, u, 0) such that {p — p, u, 0) belongs 
to 

~ K^L K I 1 iv 1-, ~ I 3 iv ,9 



and satisfies 



||(p -p, u,0)||e ^ M(||po -p|Ljv_5 iv + ||uo|Ljv_3,a:_i) 



^2.1 ^2.1 



X where M is independent of the initial data and the hybrid space -82^ 

^2,\ n ^20 /or Si^S2. 

The paper is organized as follows. We recall some Littlewood-Paley theo- 
ries for homogeneous Besov spaces and give the definitions and some properties 
of hybrid Besov spaces in the second section. In Sections 3-4, we are dedicated 
into reformulation of the system and proving a priori estimates for a linearized 
system with convection terms. In Section 5, we prove the global existence and 
uniqueness of the solution. 



2 Littlewood-Paley decomposition and Besov spaces 



Let ip : 
such that 



[0, 1] be a radial smooth cut-off function valued in [0, 1] 



1, lei ^ 3/4, 

smooth, 3/4 < |^| < 4/3, 

0, lei > 4/3. 
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Let ip{^) be the function 



¥.(0:=^(e/2)-^(0- 

Thus, ip is supported in the ball e : |^| ^ 4/3}, and (f is also a smooth 
cut-off function valued in [0, 1] and supported in the annulus G : 3/4 ^ 
1^1 ^ 8/3}. By construction, we have 

fcgZ 

One can define the dyadic blocks as follows. For k E Z, let 

A,/ := ^-V(2-'=0^/- 

The formal decomposition 

/ = EA./ (2.1) 

is called homogeneous Littlewood-Paley decomposition. Actually, this decom- 
position works for just about any locally integrable function which has some 
decay at infinity, and one usually has all the convergence properties of the 
summation that one needs. Thus, the r.h.s. of (2.1) does not necessarily con- 
verge in t5^'(M^). Even if it does, the equality is not always true in S^'{M.^). 
For instance, if / = 1, then all the projections A^f vanish. Nevertheless, (2.1) 
is true modulo polynomials, in other words (cf. [5,17]), if / G t5^'(]R^), then 
Ek&^kf converges modulo ^[M^] and (2.1) holds in ^'(M^)/^[M^]. 

Definition 2.1 Let s e R, 1 ^ p, q ^ oo. For f e ^'(M^), we write 

II/IIb|, = E2'1|a,/|U2. 



A difficulty comes from the choice of homogeneous spaces at this point. 
Indeed, || ■ ||^|^ cannot be a norm on {/ G ,y{R^) : ||/||^|^ < because 
— means that / is a polynomial. This enforces us to adopt the 
following definition for homogeneous Besov spaces (cf. [4]). 

Definition 2.2 Let s G M and m = — -|- 1 — s] . // m < 0, then we define 
5| i(M^) as 

Ki = {/ e ^'(K^) : II/IIbi andu = Y./\uf m ^'(M^)}. 
^ ' fcez ^ 

Ifm'^O, we denote by £Prn the set of two variables polynomials of degree less 
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than or equal to m and define 

Bi, = {/ e y'{R'')/^m : Wfhs ^<oo and u=J2 Afc/ m ^'(R^)/^,„|. 

For the composition of functions, we have the following estimates. 
Lemma 2.3 ( [4, Lemma 2.7]) Let s > andu E B^^^nL"^. Then, it holds 

(i) Let F e wj^l~^^'°°(R^) with F(0) = 0. Then F{u) e By. Moreover, 
there exists a function of one variable Co depending only on s and F, and such 
that 

\\Fiu)\\B^^^<Coi\\u\\L<^)\\uUs^^. 



(ii) Ifu,ve B,]„ {v-u)E Bl, fors G (-f , f ] and G e (K^^) 
satisfies G'{0) = 0, then G{v) — G{u) E B21 and there exists a function of two 
variables C depending only on s, N and G, and such that 

\\G{v) -G{u)\\^s^ ^ C{\\u\\l^,\\v\\l^) [\\u\\ n + \\v\\ N \\v-u\\j^s^. 



We also need hybrid Besov spaces for which regularity assumptions are 
different in low frequencies and high frequencies [4] . We are going to recall the 
definition of these new spaces and some of their main properties. 

Definition 2.4 Let s, t eR. We define 

Let m = — [y + 1 — s], we then define 

^2;i(K^) = {/ e ^'(M^) : ll/b-* < 00} , zfm < 0, 
S0(R^) = |/ e ^'(M^)/^^ : ll/II^M < 00} , zfm ^ 0. 
Lemma 2.5 We have the following inclusions for hybrid Besov spaces. 

(i) We have B^fi = B^^. 

(ii) Ifs^t then B^'l = B^^ n i3\-y. Otherwise, B^'l = 5| ^ + B^^. 

(iii) The space ^g'l coincides with the usual inhomogeneous Besov space 

^2,1- 
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(iv) Ifsi ^ S2 and U ^ ts, then B^'i^' ^a'i*'. 

Let us now recall some useful estimates for the product in hybrid Besov 
spaces. 

Lemma 2.6 ( [4, Proposition 2.10]) Letsi, Sa > and f, g E L^^nB^y^ 
Then fg G B'^^f' and 

Let Si, S2,ti,t2 ^ Y ^'^^^ thatmm.{si+S2,ti+t2) > 0, / G -B2,\'*^ andg G B2^f^. 
Then fg G ^s',^^^^-''*^^*^-' an(i 



For a, /3 G M, let us define the following characteristic function on Z: 



[A ifr^l. 
Then, we can recall the following lemma. 

Lemma 2.7 ( [4, Lemma 6.2]) Let F he an homogeneous smooth function 
of degree m. Suppose that —N/2 < si,ti,S2,t2 ^ 1 + N/2. The following two 
estimates hold: 

|(F(D)A,(vVa),F(D)Afea)| 

< c,2-'=(^^^'^^W--)||v||^,Ja||.n.2||F(D)A,ay^ 

-^2,1 

|(F(D)A,(v ■ Va), Afc6) + (Afc(v ■ V6), F(Z})Afca)| 

+ 2-S'^^^^^^W-™)||abn.2||A,6|U2), 
where (-, ■) denotes the L'^-inner product, the operator F{D) is defined by 
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3 Reformulation of the Original System 



Let p = p — p. Then (1.1) can be rewritten as 

Pt + u ■ Vp + p divu = —p divu, 
U( + u ■ Vu — p^^pAu — p^^(/i + A) V divu + Vp — 



- ,/ . (/xAu + (/i + A)Vdivu), ^^'^^ 
p{p + P) 



_ = p. 



Denote A'z := ^"^1^1*^2 for all s G M. Let c = A"Mivu be the 
"compressible part" of the velocity and I = curlu be the "incompressible 
part". Then, we have 

u = -A- Vc - A-^ divl, 
since divdivl = 0. In fact, 

divl = A~Miv curlu = A"^(Au - Vdivu), 

which yields 

divdivl = A~^ div(Au - V divu) = A~^( divAu - A divu) = 0. 
Moreover, 

curl divl = AI. 

The first equation in (3.1) is changed into 

Pt + u ■ Vp + pAc = — p divu. (3.2) 

For the 2nd equation in (3.1), applying A~^ div and A~^ curl to both sides, 
respectively, we get 



' ct - p-\2p + A)Ac - Ap - A-^ 

^^{pAu + (/i + A)V divu) j. , 

(3.3) 



-A~^ div J u ■ Vu + — A— (/iAu + (p + A)V divu) I 
I pip + P) J 



It - p- Vai 

= -A~^ curl |u ■ Vu + -^-^—^{pAu + (p + A) V divu) | , 

where we have used the fact 

curlV/ = {d,dj - d,d,f),, = = 
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for any function /. 

Because the first equation of (3.3) involves A^^p, we denote h = A~^p. 
Then, we have 



( ht + A'\u ■ V Ah) + pc = F, 

cj + u ■ Vc - p^\2p + A)Ac - A^h-h = G, 

It - p"VAI = H, 
^ u = -A'^Vc - A"^ divl, 



(3.4) 



where 



F = -A-^{Ah divu) , G = u ■ Vc - A"^ div J, H = -A'^ curlJ, 



J = u ■ Vu 



P 

p{p + P) 



(pAu + (p + A) V divu). 



The third equation is, up to nonhnear terms, a mere heat equation on I. 
We therefore expect to get appropriate estimates for the incompressible part 
of the velocity via the following lemma. 

Lemma 3.1 Let s G M, r G [1, +C)o], and u solve 

iut- p^^pAu = /, 
\ m(0) = Uq. 

Then there exists C > depending only on N, p~^fi and r such that, for all 
< T < +00, 

Moreover, u G C([0, T]; 5^^). 

For the first two equations, which is a linear coupling system, we can use 
the following lemma. 

Proposition 3.2 Let [h, c) be a solution of 
r ht + A-\^^■WAh)+pc = F, 

\q + v Vc-p"i(2/i + A)Ac-A2/i-/i = G, ^ ' ' 

on [0,T), (3 - N)/2 < s ^ {N + l)/2 and V{t) = Jo \\^{t)\\ ^N^.dr. The 
following estimate holds on [0, T) : 



' B 2 

^2,1 



\\h{t)\\ s + ||c(t)|| _i + [\\\hir)\\ s + ||c(r)|| +3)dr 

■°2,1 ^2,1 "'0 ■°2,1 ■°2,1 
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^Ce^^w(||M0)||_,.,3 + ||c(0)|| 

+ Te-^^W (\\F{r)\\ . + ||G(r)l| dr 

JO y i32,l ^2,1 / 

where C depends only on N and s. 

Let us define the functional space for 2 — N/2 < s ^ N/2 + 1: 

||(/.,u)|U. = + (3.6) 



When the time variable t describes a finite length interval [0, T], we will denote 
by and || ■ ll^;^ the corresponding spaces and norms. 



4 The estimates for the linear model 



This section is devoted to the proof of Proposition 3.2. Let {h, c) be a 
solution of (3.5) and denote / := e~^^^*^/ for any function /. Then the 
system (3.5) can be transformed into the following form: 



ht + A'\^r■VAh) + pc = E-KV'{t)h, 

5t + V ■ V£ - p"^(2^ + A)A£ - A^/i - /i = G - A'y'(t)c. 



Applying the operator to the system (4.1) and denoting := A^f, 
we have the following system 



dth + A"'(v ■ VAh)k + pck = Fk- KV'{t)hk, 

dtCk + (v ■ Vc)k - p-\2fx + \)Ack - K^hk -hk = Gk- KV'{t)ck, 
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4-1 Low frequencies {k ^ 0) 



Taking the scalar product of the first equation of (4.2) with hk, of the 
second equation with c^, we get the following two identities: 



( I d ..r 
2di 

1 d 



2dt 



\\hk\\h + P{^k, hk) = {Fk, h) - KV'{t)\\h4l, 

-(A"i(vVA/i)fc,M, 
W^kWh + p~\2i2 + X)\\Ack\\l2 - (A/ifc,A5fc) - ihk,Ck) 
= (Gfc, - KV'{t)rc4l, - ((v ■ V£)fe, Ck) 



(4.3) 



Now we want to get an equality involving Ahk. To achieve it, we take scalar 
product of the first equation of (4.2) with A^/i^, A^/i^ and A^c^ and of the 
second equation with A'^hk and then sum the last two resulting equalities, 
which yields, with the Plancherel theorem, that 



r 1 d 

2dt 



dt 



\\Ah\\l2 + p{Ack, Ahk) 

= {AFk,Ahk) - KV'{t)\\ AhWh - ((v ■ VA/i)fc, A/ifc), 

(A2/ifc,5fc) + p||A£fe||2, +p-i(2^ + A)(A2/ifc,A2efc) - \\A^hk\ 

- ||A/ifc||i. = {AFk.Ack)-2KV'{t){A^hk,Ck) 

- ((v ■ VA/i)fc, A5fc) + (AGfc, A/i,.) - (A(v ■ V5)fc, A/i^) 



IP- 



(4.4) 



A linear combination of (4.3) and (4.4) yields 



1 d 
Ydt 



-ll^fclli^ + -||A/ifc||i2 + ||£fc||i2 - 2Ki{A'^hk, Ck) 



P 



P 



+ K.WA^hkWh 



+ K,\\Ahk\\h + [p~\2p + A) - pK,]\\A~Ck\\i2 - p~\2fi + X)K^{A'hk, Ack) 
{Fk, hk) - KV'{t)\\hk\\l2 - (A-i(v ■ VAh)k, hk) 
1 



P 



{AFk,Ahk) - KV'{t)\\Ahk\\i2 - ((v ■ VAh)k,Ahk) 



+ {Gk, Ck) - KV'{t)\\ck\\l2 - ((v ■ Vc)k, Ck] 



{AFk, A~Ck) + {AGk, Ahk) - 2KV'{t){A%, Ck) 
- ((v ■ VAh)k,Ack) - (A(v ■ Vc)k,Ahk) 



Noticing that ||A/ifc||i2 ^ |2''||/ifc||i2 ^ |||/2,fc||^2 for A; ^ 0, we have 



(A^/ifc, Cfc) 



32Mi„ ~ „2 1 



17. I|2 
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32M2 iia2?'ii2 1 iiA~ii2 
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2M, 



Thus, we have to choose Ki, Mi and M2 satisfying 



2/i + A - p^Ki - — — 



> 0. 



Hence, we can take 



Ml 



4Vp 



5p . f p(2/i + A) 1 

iii = mm 



16(2/i + A) 



^p3 + 2(2/i + A)2' 8v/r 



Denote for k ^ 



P P 



Then, there exist constants C3 and C4 such that 



Thus, there exists a constant c such that for k ^ 



1 c/ 



+ ic2^' + KV')ai 



2dt 
1 



P 



{FkM) - {A-\-i ■VKh)kM) 
1 



P 



{AFk,Ahk) - ((v- VA/i)fe,A/i, 



+ {Gk, Ck) - ((v ■ Ve)fc, - Ki 



{AFk,A~Ck) + {AGk,Ahk) 



- ((v ■ VAh)k, Ack) - (A(v ■ Ve)fc, A/iA 
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4-2 High frequencies {k > 0). 



Taking the scalar product of the first equation of (4.2) with Ahk, of 
the second equation with Ac^, we get the following two identities: 



= {Fk,Ah) - KV'{t)\\A-2hk\\h - ((v ■ VAh)k, hk), 
~\\Ahk\\l2 + p~\2i2 + \)\\Ahk\\l2 - iA^hk,Ack) - iAhk,Ck) 
= (A^Gfc, A^£,) - KV'{t)\\A^~Ck\\h - ((v ■ Ve)fc, A£fc). 



(4.6) 



Now we want to get an equality involving A^hk. To achieve it, we take 
scalar product of the first equation of (4.2) with A^hk, A^hk and A^Ck and of 
the second equation with A^hk and then sum the last two resulting equalities, 
which yields, with the Plancherel theorem, that 



r 1 d 

2dt 



1 d 
2di 



d 
dt 



\\A^hk\\l2 + p{Ack,A^hk) 
= {Ah.A^hk) - KV\t)\\A^hk\\l2 - ((v ■ VA/i)fc, A^/i,), 
\\Aihk\\l2 + p{A^~Ck, A'^hk) 
= {A'Fk,A%) - KV'mAlhkWl^ - (A(v ■ VAh)k,A'hk), 

(A^hk^Ck) + p\\Ahk\\l2 + p~\2i^ + X){A^hk,A^ik) 

-||At/i,,||i.-||Ai/i,||i. 
= (A^Ffc, A5,) - 2KV'{t){A%, c^) - ((v ■ VAh)^, A^Ck) 
+ {AGk,A^h) - {A{^^ ■Vc)k,A%). 



(4.7) 



A linear combination of (4.6) and (4.7) yields 



1 d 
2di 



-\\A-2hk\\i2 + -\\A-^hk\\i2 + p-'{2fi + X)K2\\A^^hk\\i2 



P 



P 



+ \\Ahk\\l2-2K2{A^hk,Ahk) 
[p~\2fx + A) - pK2]\\ Al~c4l2 + K^WA'^hWh + ^2||At/ifc||i. 
iFk,Ah) - KV'{t)\\A-2h\\l2 - ((v ■ VAh)k, hk] 
1 

+ - 
PL' 

+ p-^{2p + X)K2 



{AFk,A%) - KV'it)\\A^^h\\l2 - iiy^-VAh)k,A%) 



{A'Fk,A'hk) - KV'{t)\\A-2hk\\l2 - (A(v ■ VAh)k,A^h, 
(A^Gfc, Ahfc) - KV'{t)\\A-2c4l2 - (A5(v ■ V5)fc, A^g^) 
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Noticing that 



- (A(v ■ VAh)k, Ack) - (A(v ■ Ve)fc, A%) 



(A2/i,,A5cfc) ^ —\\A-^hk\\l2 + ^IIA^Cfclli^, 



we have to choose K2, M3 such that 

< < M3 < p-2(2/i + A). 
For example, we can take 

J 0-2 ' ^ 



4p' 



7^2 



Denote for k > 

2 1 
at. :=- 



1 



||A2/i,||i, + ^WA-^hWh + p"^(2/i + A)ir2||A2/i,||i, 



P P 
+ ||Ahfc||i2 -2A'2(At/ifc,Ahfc). 

Then, there are constants Ci and C2 such that 

cia^. ^ WA-^hkWh + WA'^hkWh + \\A~^~hk\\h + \\Ahk\\l2 ^ C2a^ 

Thus, by Bernstein's inequality ||A^/ifc||/^2 ^ |2'^|| A/;,fe||i2, there exists a con- 
stant c such that 

1 d 



2dt 
1 

1 

+ - 
P 



al + {c + KV')al 



{Fk,Ahk)-{{Y-VAh)k,hk) 



(AFfc,A2/^fc)-((v VA/i)fc,A% 



(A^Ffc, At/ifc) - (A2(v ■ VA/i)fc, A2/i, 



(4.8) 



+ (A2G'fc,A2£fc) - (A2(V ■ V£)fc,A2£;t) 



(A5Ffc, A2 5fc) + (A^Gfe, AS/ifc) - (A(v ■ \/Ah)k,Ack) 



(A(v- V£)fc,A2/i, 



Now, we combine (4.5) and (4.8). At this stage, we use Lemma 2.7 to 
estimate the terms involving a convection in (4.5) and (4.8), and eventually 
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get the existence of a sequence {■~fk)kez such that J2kezlk ^ 1 and 
1 d 



ai + {cmm{2'\l) + KV')ai 

||(F,G)|L._,.,3 



2dt 



5- +^'ll(/^,5)ll^.-M+i^^.-M-^ 

-°2,1 X-°2,l -"2,1 ^^2,1 



(4.9) 



where c = minfc, c). 



We are going to show that inequahty (4.9) provides us with a decay for 
h and c. We actually have a parabolic decay for c. 



4-3 The damping effect for h 



Let (5 > be a small parameter (which will tend to 0) and denote f3l 
al + 6"^. From (4.9) and dividing by Pk, we get 



d 



(3k + {cmm{2^'',l) + KV')(3; 



dt 

^C-fk2~''^"~^^ 



||(F,G)|L._,,.,3 



D 2 y D 2 



-°2,1 ^^"2,1 



+ y'||(/.,5)|| 3 

^2,1 

+ 6{cmm{2^\l) + KV). 

Integrating over [0, t] and making 6 tend to 0, we have 

i-t 

afc(t) + cmin(2^'', 1) / akir)dT 
Jo 

^afc(0) + C2-^'(-i) r7,(r)||(F,G)|U._, 



(4.10) 



_3 



,s + * -s-l,s- 



.dr 



(4.11) 



+ / V'{t) 
Jo 



B^ ^xB^ 2 



i.s-i - Kukir) 



dr. 



By the definition of a^, we have for any A; G Z 



2*^(^-i)a2 ^ 2'=(^-i) max(l, 2i'=)||^fc||i2 + 2'=(*-i) max(l, 2^\\ck\\l2. (4.12) 
Thus, we have in taking K large enough such that 



E 



C2-'^^-^hk{r)\\{kc)\\ 



.-i,.+ 3 -.s-i,s-i - i^afc(r) 

^2,1 ^^2,1 



^ 0. 
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Multiplying both sides of (4.11) by 2^^^ ^\ According to the last inequality, 
and due to (4.11) and (4.12), we conclude after summation on k in Z, that 



IIM^)IL.-i..+a + ||£(t)IL.-M 1 + c / ||/i(r)|L,+,„+3 



^2,1 -°2,1 -"2,1 

+ ^ rc2^-(^-^)min(2f,l)||£,(r)|U2cir 

° (4.13) 



^ll(MO),c(0))|| 3 ^_,,,_i 



+ ||(F,G)(r)|| 3 _i._irfr. 



^.^ T/ie smoothing effect for c 

Once stated the damping effect for /i, it is easy to get the smoothing effect 
on c. Since (4.13) implies the desired estimate for low frequencies, it suffices 
to prove it for high frequencies only. We therefore suppose in this part that 
k>0. 

Define 6k = \\ck\\L^. By the previous inequalities and using Lemma 2.7, 
the second equality of (4.3) yields, for a constant c > 0, that 

+ ^fcl^'(t)(C7fe2-'=(^-i)min(l,2-t)||5||_^_^^^_, _A'^,). 

^2,1 

Using j3l = 6l + (5^, integrating over [0, t] and then having 6 tend to 0, we infer 

ek{t) + c2^''ek m{0) + C f'wGkh^dr + C f\^''\\hk{r)\\L^dT 

Jo Jo 

+ C f V'(r)7fc(r)2-^(-i) min(l, 2"i)||£(r)|| ^._i,._i cir. 

Jo ^2,1 

Therefore, we get 

^ 2^^(-i) II c,(t) IU2 + c r ^ 2^'(^+i) II 5fc (r) |U2rfr 

fe>0 ° fc>0 

^l|£(0)IL.-i.-^ +C f \\G{t)\\ .dr + C fY. 2'=(^+t)||/.,(r)|U2rfr 

+ Cy(t) sup ||£(t)|| 

[o,t] ^2,1 
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Using (4.13), we eventually conclude that 




^{C + CV{t))(\\hm . + ||5(0)|| i) 

\ ^2,1 ^2.1 

+ A||%)|L.-M+i + l|G'(r)|| ,._iMr). 

JO i'2.1 ^2,1 ' 

Combining the last inequality with (4.13), we complete the proof of Proposi- 
tion 3.2 as long as we change the functions {h, c, F, G) back into the original 
ones (/i, c, F, G). 



5 A global existence and uniqueness result 



This section is devoted to the proof of Theorem 1.1. The principle of the 
proof is a very classical one. We shall use the classical Friedrichs' regularization 
method, which was used in [2,3,16] for examples, to construct the approximate 
solutions (/i",u") of (3.4). 



5.1 Building of the sequence (/i",u'^)„gN 



Let us define the sequence of operators {^n)n&n by 
/J- := ^-il^(i,„)(0=^/. 
We consider the approximate system: 

' Kl + /•„A-i(/-„u" ■ VA^„/i") + p/-„c" = F\ 
< + ■ V/-„c") - -p~\2^ + A) A/-„c" 

- A2/-„/i" - = G 

u" = -A- Vc" - A-i divF, 
^ (/i",c",r)(0) = (/i„,A-Mivu„,A-icurlu„), 

where 

= -^nA-\A /:„/i"div /'nU"), 

= M^nu^ ■ V Ac") - divJ", 
if" = -^A-^curlJ", 
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pC(A/-„/i" + p) 
where C is a smooth function satisfying 



(/iA^„u" + (/i + A) V div/-„u'^ 



p/4, |s| ^ p/4, 
s, p/2 ^ |s| ^ 3p/2, 

7p/4, |s| ^ 7p/4. 
smooth, otherwise. 



We want to show that (5.1) is only an ordinary differential equation in 
X X L^. We can observe easily that all the source term in (5.1) turn 
out to be continuous in x x L^. For example, we consider the term 
^nA~ div fe"+p) • Plancherel's theorem, Hausdorff- Young's inequal- 
ity and Holder's inequality, we have 



11^"^ ^" C(AA/^- + p) 



L2 = l|li.(ln)ler'e-^ 

\ ?7, ' ' 



A/-„/z"A/-^u" 
C(A^„/i" + p) 



i2 



C(A/-„/i« + p) 



4 4??^ 
^^||A^„/^"|Uo.||A^„u"iU. ^ —II |e|lB(in)^/^"1Ui||u"|U2 
p p ' 

4n'5'"'"^ 

^ ILnll ||„n|| 

~ — = — IIl2||u 11^2. 

P 

Thus, the usual Cauchy-Lipschitz theorem implies the existence of a strictly 
positive maximal time T„ such that a unique solution exists which is continu- 
ous in time with value in x x L^. However, as = we claim that 
^n{h^^ c",I") is also a solution, so uniqueness implies that ^n{h^, c",I") = 
(/i", c", !"■). So (/i", c", !"■) is also a solution of the following system: 

+ /■„A-^(u" ■ VA/i") + pc" = Fi", 
cr + ■ Vc") - p~'(2/i + A) Ac" - A'/i" - /i" = 

IJ^ - p-VAF = iff, (5.2) 
u" = -A- Vc" - A-MivF, 
[ (/i",c",r)(0) = (/i„,A"Mivu„,A"icurlu„), 



where 



Fi" = -/-„A-i(A/i'^divu'^), 

= /•„(u" ■ Vc") - Jf^k-' div J", 
ifr = -/-„A"icurlJ", 
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Jr = ^„(u" ■ Vu") + _^^^^„^_^ (Mu" + (/i + A)Vdivu"). 
The system (5.2) appears to be an ordinary differential equation in the space 




Due to the Cauchy-Lipschitz theorem again, a unique maximal solution exists 
on an interval [0,T*) which is continuous in time with value in L^^x L\x L"^. 



5.2 Uniform hounds 



In this part, we prove uniform estimates independent of T < in i?" 
for (/z", u"). We shall show that T* = +00 by the Cauchy-Lipschitz theorem. 
Define 

E{0) :=||A-^(po - 4.3,4^, + lluoll 4_3,4_i, 

^2,1 ^2,1 



E{h,u,t) :=||(/i,u) 



E_ 



N . 



fn := sup [t e [0, T*) : u", t) ^ ACE{0)} , 

where C corresponds to the constant in Proposition 3.2 and A > max(2, C~^) 
is a constant. Thus, by the continuity we have T„ > 0. 

We are going to prove that T„ = for all n G N and we will conclude 
that T* = +00 for any n G N. 

According to Proposition 3.2 and Lemma 3.1, and to the definition of 
{hn,Un), the following inequality holds 



C||u"|| jv+i . 



T ^ 2,1 

Un JV 3 JV 1 + LfT JV 3 JV 

- Vc" JV„3JV_, + JV _ 3 JV 



Therefore, it is only a matter of proving appropriate estimates for F", 
J" and the convection term. The estimate of is straightforward. From 
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Lemma 2.6, we have 



IIK^II iV 3 N,, =||A/l'^diVu"|| N 5 N 



^C||A/i"|| iv_o iv II divu"|| N^i N 



^C\\h''\\ N 3 N., ||U" 

<CE2(/i",u",T). 



JV , 1 AT 



2 1 -in ,,n 



With the help of Lemma 2.6 and interpolation arguments, we have 

U ■ VC iV_3 iV_-, 

^C||u"|| iV 1 ^ 1 ||VC"|| iV 1 N 



U JV 3 iV -I U JV , 1 JV 



In the same way, we can get 



|u"-Vu"|| iv 3 JV , ^CE\h^,u'\T). 



To estimate other terms of J", we make the following assumption on 

2CiACE{0) ^ p, 

■ K 

where Ci is the continuity modulus of the embedding relation -82^1 ( 
L°°(R^). If T < fn, it implies 



l|A/l"|Uoo([o,T]xEiV) ^Ci||/l"i|^^^^^^^+i^ ^ Cl||/l"||^^^^^^^.3,^+,^ 

^C,ACE{0) ^ ip. 



which yields 



A/i" + p G [-P, -p] and C(A/i" + p) = A/i" + p. 
2 2 
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From Lemma 2.6 and Lemma 2.3, we obtain 
— (/iAu" + + A)Vdivu" 



II ^ , , ^ V A^"" I VP- I ^i/vu-ivu ;|| iV_3]V_-, 

" p(A/l" + p) i'M^a'S ^'^ ) 
<C\\— II JV U N ,1 AT ,-, 

^C||A/l"|| AT U AT I 1 JV , , 

^c^ll/i"ll 

JV 3 JV ,-, U iV I 1 JV 11 

^CE\h'',u'',T). 
Thus, we get 

||(/i",u")|| JV ^(7e^^'^(°Hl + C^'C'2E(0)]E(0). 

So we can choose -E'(O) so small that 

1 + CA^C^E{Q) ^ e^'^^^o) < 2C^ACE{^) ^ p, (5.9) 

which yields ^iU" , u")||£;i ^ ^AC'E(O) for any T < f;. It follows that T„, = 

T„*. In fact, if f„ < T„*, we have seen that u'^, T^) ^ ^AC'E(O). So by 

continuity, for a sufficiently small constant a > we can obtain u", T„ + 

cr) ^ y4C*-E'(0). This yields a contradiction with the definition of T„. 

Now, if f„ = T„* < cx), we have obtained u", T„*) ^ A(5e(0). As 

< and ||u"||^^^(gO,.) < oo, it implies that ||/i"||l^.(l2) < oo 

and ||u"||j;^^,(j;^2) < oo. Thus, we may continue the solution beyond T^* by the 
Cauchy-Lipschitz theorem. This contradicts the definition od T^. Therefore, 
the approximate solution (/i", u"')„gN is global in time. 



(5. 



5.3 Existence of a solution 



In this part, we shall show that, up to an extraction, the sequence u"')^^^!} 
converges in ^'(M"*" xM^) to a solution {h, u) of (3.4) which has the desired reg- 
ularity properties. The proof lies on compactness arguments. To start with, we 
show that the time first derivative of (/i", u") is uniformly bounded in appro- 
priate spaces. This enables us to apply Ascoli's theorem and get the existence 
of a limit {h, u) for a subsequence. Now, the uniform bounds of the previous 
part provides us with additional regularity and convergence properties so that 
we may pass to the limit in the system. 

It is convenient to split (/i", u") into the solution of a linear system with 
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initial data (/i„,u„), and the discrepancy to that solution. More precisely, we 
denote by (/i^, u^) the solution to the linear system 

' dthl + pA"^ divu^ = 0, 

< dtul - p-^fiAul - p-\fi + A)Vdivu" + VAhl + VA'^hl = 0, (5.10) 

and (/i", u") = {h"" - hi, u" - u^). 

Obviously, the definition of u„) entails 

_ iV_ 3 iV_|_-|^ ^ N__3 N__l 

hn — > A" (po — p) in i?2^i ^' ^ , u„ — i> Uq in B2^i ^' ^ as n — oo. 
The Lemma 3.1 and Proposition 3.2 insure us that 

(/i2,0^(/iL,Ui)ini?^, (5.11) 
where [hi, vll) is the solution of the linear system 
dthl + pA""*^ divui = 0, 

9tUi -p~VAui -p"^(/" + ^)Vdivu + VA/ii + VA^i/iL = 0, (5.12) 
{hL, UL)t=0 = (A"^(po - p), Uo). 

Now, we have to prove the convergence of {h'^,u"'). This is of course a 
trifie more difficult and requires compactness results. Let us first state the 
following lemma. 

Lemma 5.1 ((/i", u"))„gN is uniformly bounded in 

1 • 1 • i^_3 

Proof. Throughout the proof, we will note u.b. for uniformly bounded. We 

. JV_3 

first prove that dth^ is u.b. in (L^ + L°°)(M"'', B21 ^ ), which yields the desired 
result for h. Let us observe that /i" verifies the following equation 

= _ ^„A-i(A/i"divu") - /•„A-i(u" ■ VA/i") 
- pA~^ divu" + pA~^ divu^. 

According to the previous part, (/i"')„gN is u.b. in L°°{B2^i ^) and (u")„gN is 

. K 

U.b. in L^(i?2^i) in view of interpolation arguments. Thus, ^„A~^(A/;," divu"), 
^„A-i(u" ■ VA/i'^), pA-Mivu'^ is u.b. in L'^{B^^i '). The definition of 
obviously provides us with uniform bounds for A~^ divu^ in L°°(i?2^i ^), so 
we can conclude that dth^ is u.b. in (L^ + L°°){B2^^ ^). 
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Denote = A-MivuJ, c"" = A-Mivu", II = A'^ curluj and P = 

A"^ curlu'^. Let us prove now that dtc"- is u.b. in (Ls + i?2^i ^) and 

that Stl*^ is u.b. in Ls (M+; i^g^-,^ which give the required resuh for u" by 
using the relation u*^ = — A~^ Vc" — A~^ divP. 

Let us recall that 

a^c" =p-\2fi + A)A(c" - cl) + A2(/i" - hi) 

+ (/i" - - ^„A-MivJ", 
dtr =p- VA(P - 12) - A^"^ curlJ'^. 

Results of the previous part and an interpolation argument yield uniform 

bounds for u" and c" in ^3(^2^^^') n L^iB^^). Since /i" is u.b. in L'^iBii ) 
4 • —+- 

and is u.b. in Ls (i?2\ ^ ), we easily verify that A(c" — c^) and ^n-^~^ divJ" 
are u.b. in L3(i?2^^ ^). Because /i" is u.b. in L°°{B2^i A^/i" is u.b. in 
L°^{Bi^i We also have K^hl u.b. in L'^iBi^i In addition, /i" and /i^ are 
U.b. in L°°(52^i gQ gj^aiiy u.b. in (L3 + L°°)(M+; ^2^1 r^j^g 

case of goes along the same lines. As the terms corresponding to (/i" — h^) 

4 . iY._3 

do not appear, we simply get c^iP u.b. in L3(i?2^i ^). □ 

Now, we can turn to the proof of the existence of a solution and use 
Ascoli theorem to get strong convergence. We need to localize the spatial space 
because we have some results of compactness for the local Sobolev spaces. 
Let (Xp)pGN be a sequence of C^(M^) cut-off functions supported in the ball 
B{0,p + 1) of and equal to 1 in a neighborhood of B{0,p). 

For any p E N, Lemma 5.1 tells us that ((XpP") XpU"'))nGN is uniformly 

• iy._3 

equicontinuous in C(]R"*"; (i?2^i ^)^^^). 

Let us observe that the application / i— > Xpf is compact from i?2^^ ^' ^ 

. iL — l ~ £L — 1 K — i ■ K — S. 

into -82^1 ^ , and from 52^i ^' ^ into i?2^i ^- After we apply Ascoli's theorem 
to the family {{Xph"', XpU"))neN on the time interval [0,p], we use Cantor's di- 
agonal process. This finally provides us with a distribution (h, u) belonging to 

C(M"'"; (-82^1 2)^+^) and a subsequence (which we still denote by ((p", u'^)^^^) 
such that, for all p eN, we have 

{Xph^,Xp^n ^ iXpKXp^) as n ^ +00, in C([0,p]; (B^^, ^)^+^). (5.13) 
This obviously infers that (/i", u") tends to (h.u) in ^'(M+ x M^). 

Coming back to the uniform estimates of the previous part, we moreover 
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j;et that {h, u) belongs to 

_JV_3JV,-, ~ N _3 N \ / _iV|lJV|-| 



and to Ci/2(^+.5^2^ 2) X (cV4(K+;52^^ 2))7V^ 

Let us now prove that (/i, u) := (/i^^jU^) + (/i, u) solves (3.4). We first 
observe that, according to (5.1), 



' + ^^"'(u" ■ VA/i") + pc" = -^„A-i(A/i"divu"), 

(5.14) 



< + /■„(u" ■ Vu") - p^VAu" - p^\fi + A)Vdivu" + AV/i" 



A/?" 

A-iy/,- = - A ^ _^ (/.Au" + (/. + A)Vdivu") 



The only problem is to pass to the limit in ^'(M+ X M^) in nonlinear 
terms. This can be done by using the convergence results stemming from the 
uniform estimates and the convergence results (5.11) and (5.12). 

As it is just a matter of doing tedious verifications, we show, as an 
example, the case of the term p[Ah"+p) ■ Denote L{z) = z/{z + p). Let 
e e Co=^(M+ X M^) and p G N be such that supp^ C [0,p] x 5(0, p). We 
consider the decomposition 

^ ^A/i"Au" ^A/iAu 



^ p(A/i" + p) p{kh + p) 
= -p~^f^[e{l - L(A/i"))xpA/i-XpA(u2 - u,.) 
+ ^(1 - L(A/i"))xpA/i"xpA(xp(u" - u)) 

+ ^(1 - L(A/z"))(xpA(/i" - h))An - M/iXpAu(L(xpA/.") - L(xpA/i))] 
, ^ eAhAu 

The last term tends to zero as ^ +oo due to the property of As 
eL{Ah'') and A/i" are u.b. in L°°{Bi^^) and tends to in Ls{Bl^ 2), the 
first term tends to in Ls{B2i ^). According to (5.12), Xp(u" ~ u) tends to 

4 . iV_3 ' _^ . iV_3 

zero in L3[B2i ^) so that the second term tends to in Li{B2i ^). Clearly, 

' . JV ' ■ E. 

XpA/i" ^ XpA/i in L°°{Bi^i) and L(xpA/i") L{xpAh) in L°°(L°° n ^2^1), 

4 . iV._3 

SO that the third and the last terms also tend to in Li{B2i ^). The other 
nonlinear terms can be treated in the same way. 



~ K — l — +1 



We still have to prove that h is continuous in i?2^i ^ ' ^ and that u 
belongs to C(M"''; B21 ^' ^ ^). The continuity of u is straightforward. Indeed, 
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u satisfies 



Ut = — u ■ Vu + p /_iAu + p (/i + A) V divu 

- AV/i - A~ V/i - — -(/iAu + (/i + A) V divu) 

p[kh + p) 



_ JV 3 JV 



and the r.h.s. belongs to (L^ + L°°)(M"''; B2^i ^' ^ ). We have already got that 
h G C(M+; ^2^1 2). Indeed, /i* G L°°(R+; ^2^1 2) from the equation 

ht = -A~^ div(A/iu) - pA~^ divu. 

Thus, there remains to prove the continuity of h in i?2^i • 

Let us apply the operator dk to the first equation of (3.4) to get 

dtAkAh = -Afe(u ■ VA/i) - pA^divu - Afc(A/idivu). (5.15) 

Obviously, for fixed k the r.h.s. belongs to ^^^^^(IR"''; L^) so that each AkAh is 
continuous in time with values in L^. 

Now, we apply an energy method to (5.13) to obtain, with the help of 
Lemma 2.7, that 

~\\AkAh\\l2 ^ C\\AkAh\\L2(-fk2'''^\\Ah\\ .n\\u\\ ,n^, + ||Afedivu||i2 

jLi (Jib \ B2 2 ^2 1 

+ ||A;^,(A/idivu)||i2), 

where J2k£zlk ^ 1- Integrating in time and multiplying 2 ~, we get 

2'=(f Afc/i(t)|U. ^2'=(f -*-^)|| AfcA-i(po - p)\\l^ 

+ C f U\\Kt)\\ .,M{t)\\ 

JO \ -«2,i 

+ 2^'(f II Afcu(r) IU2 + 2^T II A,.(A/i divu) (r) rfr. 
Since h G L'^iB^y ), u G L^B^y ) and A/idivu G L^B^^^, we can get 

E 2^'(f +^)|| A.Mt)iu. <||po - Pll^^ + (1 + .OII-ll,.^f 

fcez 2,1 i.-°2,i i ^ (^2.1 I 

+ ||A/idivu|| N < 00. 

Thus, I]|fc|i;Ar converges uniformly in L°°{W^; i?2^^ ) and we can conclude 
that h G C(M+;52V )• 
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5-4 Uniqueness 



Let (/ii,ui) and (/i2,U2) be solutions of 

' ht + A"\u ■ VA/i) + pA"^ divu = -A-\Ahdwu), 
ut + u- Vu-p~V^u-p~^(/i + A)Vdivu + AV/i + A"V/i , . 

— (/iAu + (/i + A)V divu). 



p(A/i + p) 



in i?^ with the same data (A^^(po — p), uq) constructed in the previous parts 
on the time interval [0,T]. Denote {6h,6u.) = (/i2 — hi,U2 — ui). From the 
(5.16), we can get 



f dt6h + A"^(u2 ■ VA/i2) + pA~^ div5u = Fi, 
dt6u + U2 ■ V5u + p~^fiA5u — p^^ifi + A)V div5u 

+ AV6h + A-^V6h = F2 

[{6h,6u) = {0,0), 



(5.17) 



where 



Fi = - A'\6u ■ Vhi) - A~^(A5/idivu2) - A'\Ahi divSu), 

F2 = - (5u ■ Vui - _J fxASu + (/i + A) V div5u) 

p(A/i2 + p) 

+ ( A.^^- --- A/.\ - + -) (/"Aui + (/i + A)Vdivui). 
\Ah2 + p p Afti + p p) 



Similar to (3.1), we can get 
\\{5h,5M)\\ . 

^Ce ^^'''2^1 ' ll^ill ^ 3^+1 +11^211 3^ 



Noticing that 



and 



/ii, /i2 e L^{Bi^, - ^ ^ ) n L^(52^i^^' ^ ^ ) 

_]V_3 iV_-| ~ I 1 jV,-, 

Ui, U2 G L^(52^i ^' ^ ) n 4(^2^1 ^' ' ^ ) 



II^i||l°°([0,T]xR'V) ^ II^2||l°°([0,T]xR'V) ^ 
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by the construction of solutions, we have with the help of interpolation argu- 
ments 



llFlll AT 3 JV ,-, <||/il|| N 1 JV \\6u\\ N 1 JV 

+ \\6h\\ N _i JV ,-, U2 iv , 1 JV 

+ \\h\\ N^l^N \\5U\\ N^l^N 



and 



LFo JV 3 JV -I 
<||(5u|| JV 3 JV -, lluill JV.-i + (1 + ||/l2ll N.^ ] 

• ll^lll ^+1 ll^ull JV , 1 JV 

+ (||/il|| + , + ||/i2||_,,^+i l|ui| 



Thus, we obtain 



C||U2|| 



where limsup^_,Q+ Z(T) = 0. 



Supposing that 2C(l + pCi ^)y4Ci?(0) < ^ besides the conditions in (5.9), 
and taking T > small enough such that C||u2|| . jv_|_i ^ In 2 and Z[T) < i, 

we obtain 5u)|| jv = 0. Hence, (/ii,ui) = {h2,vL2) on [0,T]. 

Let Tm (supposedly finite) be the largest time such that the two solutions 
coincide on [0,Tm]. If we denote 



{Kit), u,{t)) := {Kit + T^), u,(t + T„)), z = 1, 2, 
we can use the above arguments and the fact that 

||^i||L-(R+xRiV) < and fhiL n_sn^, ^ ACEiO) 

to prove that (/ii,ui) = (/i2, U2) on the interval [0, T^] with the same as 
in the previous. Therefore, we complete the proofs. 
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